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Using the relativistic Hartree-Fock approximation, we calculate the rates of atomic ionization by 
absorption of pseudoscalar particles in the mass range from 10 to ~ 50 keV. We present numerical 
results for atoms relevant for the direct dark matter searches {e.g. Ar, Ge, I and Xe), as well as the 
analytical formula which fits numerical calculations with few per cent accuracy and may be used 
for multi-electron atoms, molecules and condensed matter systems. 
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I. INTRODUCTION 

Different lines of evidence consistently point to the ex- 
istence of cold dark matter, that comprises about 1/5 
of the current energy budget of the Universe, 5 times 
more than the usual matter (see, e.g. reviews [IHEj]). 
So far, the dark matter is seen only via its gravitational 
effects, which offers few clues of its true identity. Sta- 
ble on cosmological time scales elementary particles with 
weak or super-weak interactions to visible matter natu- 
rally arise in many extensions of the Standard Model (see 
e.g. [l|, |3[), and may comprise the bulk of dark matter. 
The last two decades have seen a many-fold expansion 
of the experimental program aimed at direct detection of 
energy deposited by the dark matter particles in the col- 
lision with atoms in low radiation environments. So far 
no positive dark matter signal has been detected, which 
places important constraints on many extensions of the 
Standard Model that predict such particles. 

Heavy stable particles (WIMPs) with masses above 
1 GeV can manifest themselves by depositing some frac- 
tion of their kinetic energy in the keV-to-100 keV range in 
the process of scattering with atoms. Many dark matter 
search strategies are tuned with such signature in mind. 
Recently, it has been pointed out that very weakly unsta- 
ble dark matter bosons in the mass range below 1 MeV 
(bosonic super- WIMPs) can be detected via its absorp- 
tion by atomic shells Q. The deposited energy is dom- 
inated by the dark matter rest mass, and leads to the 
7-like mono-energetic energy release. Several experimen- 
tal analyses exploiting such possibility have been per- 
formed @, Hi- (Note however, that the suggestion that 
the annual modulation signal observed by DAMA and 
DAMA-Libra [9( is caused by the absorption of super- 
WIMPs [lOj does not hold because inelastic processes do 
not have velocity modulation @.) Future large-scale ex- 
periments with self-shielding capabilities (llMl3l | will be 
able to improve the sensitivity to the super- WIMPs by 
several orders of magnitude. 

So far the calculations of absorption rate have been 



done using simple models of atoms which ignore the rel- 
ativistic and many-body effects [f|. This way one can 
relate the absorption cross-sections for axial, scalar and 
vector dark matter particles to the photoionization cross- 
sections. It is clear however, that atomic corrections can 
be non-negligible, and dedicated atomic calculations of 
super- WIMP absorption rates by atoms are required. 

In the present work we perform ab initio calculations of 
the ionization cross section of atoms by pseudoscalar par- 
ticles using the relativistic Hartree-Fock (RHF) method. 
For concision, we call such a pseudoscalar particle " an ax- 
ion" , stressing that the mass— coupling range considered 
in this paper does not correspond to a QCD axion. Our 
calculations are not linked to the photoionization and 
can be directly used for the analysis of the experimental 
data. We consider in particular the ionization of several 
atoms such as argon, germanium and xenon since these 
elements are used in the detectors (see, e.g. 00]). 

We have fitted the RHF results by an analytical for- 
mula. The formula works very well when axion energy is 
sufficiently large to excite an electron from the Is, 2s or 
2p states in atomic core. We checked that the discrep- 
ancy between the Hartree-Fock and analytical results is 
at the level of a few per cent for axion energy between 
10 and 50 keV and nuclear charge 18 < Z < 60. The for- 
mula can be used for any multi-electron atom and wide 
range of axion energies. The results of this work are to 
be used in the experimental analyses searching for dark 
matter using underground detectors. 



II. THEORY 

The Hamiltonian for the pseudoscalar particle a inter- 
acting with electrons can be written as [6| 



Ja 



(1) 



where energy scale parameter f a parametrizes the 
strength of the interaction; ijj is electron Dirac field. We 



2 



do not include the direct interaction of axion with pho- 
tons afl on account of strong constraints from the 7- 
ray backgrounds 6] . It is convenient to present the total 
cross section of the atom ionization by axion in a form 
where the part which needs numerical analysis is pre- 
sented as a dimensionless factor. Therefore we write 



-K{e a )al, 



(2) 



where eo is an energy scale (in our calculations eo=l a.u. 
= 27.21 eV, but it can also be any other energy unit), c 
is speed of light, v is the axion velocity in the laboratory 
frame, ao = 0.52918 x 1CP 8 cm is Borh radius, e a is axion 
energy and K(t a ) is dimensionless factor given by 



K(e a ) 



f) e £i&(2j c + l). (3) 



Here summation goes over all electron states c in the 
atomic core, 2j c + 1 is the occupation number of the 
subshell c, j c is the core electron angular momentum, 
R ce is the radial integral for the Hamiltonian ([T]) with 
two single-electron wave functions. Initial electron state 
is the state c in atomic core while final electron state is 
the state in continuum with the energy e = e a + e c . The 
integral is non-zero only between states of the same total 
angular momentum j but of the opposite parity. The 
integral can be written as 
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(/c(r)ff e (r)- ffc (r)/ e (r))dr, (4) 



where / and g are upper and lower radial components of 
the electron wave function in the spherically-symmetric 
case: 



V>(r) 



iag(r)ttji m (-n) 



(5) 



Here a is the fine structure constant (a w 1/137.036). 
The states in the continuum are normalized to the 6- 
function of the energy. Notice that we use realistic esti- 
mates of velocities of dark matter particles, v/c ~ 10~ 3 , 
which allows to good accuracy retain only the time 
derivative term in Eq. (fTJ), and the only velocity de- 
pendence is then 1/v in Eq. 

We use relativistic Hartree-Fock (RHF) method to cal- 
culate K using §S§ and ((3]). We will focus on the axion 
mass window 10 keV< e a < 50 keV, which is in the max- 
imal sensitivity range for the dark matter search exper- 
iments. Above these energies, the axion-induced Comp- 
ton like scattering ae — > ej may dominate over direct 
axio-absorption 6]. As stated before, we consider cold 
non-relativistic axions so that we can always assume that 
e a = m a c 2 , where m a is axion mass. 

In principle, formula @ can be used at any axion en- 
ergy. However, at low energy, when only outer atomic 
electrons can be excited, RHF approximation might be 



insufficiently accurate due to many-body and environ- 
mental effects not included into the RHF potential. On 
the other hand, when ionization is dominated by exci- 
tation from K and L atomic shells these effect can be 
neglected. As known, the correlation corrections to RHF 
results for inner electrons Is, 2s, 2p decrease as 1/Z where 
Z is the nuclear charge. Moreover, RHF often gives rea- 
sonable results for Is even for Z = 2. As a next step, 
the cross section can be approximated by an analytical 
formula which fits RHF results. The advantage of hav- 
ing such formulae is obvious: the results can be easily 
obtained for any atom. 

For the non-relativistic operator (but relativistic wave 
functions) the radial integral (Ql becomes (in atomic 
units) 



Rce = / (fc(r)fe(r) + a 2 g c (r)g € (r))rdr. 



(0) 



This is the same form of the radial integral which ap- 
pear in the photoionization cross section. In this case 
the expressions for the ionization cross sections by axion 
and photon are very similar and differ by a factor only. 
It was found in Ref. [6j that the cross section of atom 
ionization by non-relativistic axions can be expressed via 
the photoionization cross section cr 7 , 



<J a v 
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(7) 



Using expressions and we can present photoion- 
ization cross section in a form 



eo 



(8) 



which uses the same dimensionless function of energy 
K(e) as the axion cross section. This gives us an 
opportunity to check numerical calculations using exper- 
imental data for the photoionization. We have done this 
for the cases of krypton and xenon using the data from 
Ref. [l6| . We found that numerical and experimental 
results agree within a few per cent accuracy with an ex- 
ception of the near the threshold ionization where the 
difference is a little higher than 10%. 

Expression ([7]) helps us to find an analytical formula 
for the atomic ionization by axion absorption. Analyti- 
cal expressions for photoionization in Coulomb field can 
be found in many quantum mechanics textbooks (see, 
e -g- [13, Ell)- Using the textbook expressions for pho- 
toionization from Is, 2s and 2p Coulomb states and the 
expression (J7J) we firstly obtain the results for a single 
electron atom. Then we use additional parameters to 
fit the numerical results of the relativistic Hartree-Fock 
calculations. This way we arrive to the formulae which 
describe ionization of multi-electron atoms by absorption 
of axion: 



K, 



total 



K ls + K 2 



K> 
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(9) 
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TABLE I: Hartree-Fock energies of the core states of Na, Ar, 
Ge, I and Xe (atomic units, 1 a.u.= 27.21 eV). 
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Atom 


Na 


Ar 


Ge 


I 


Xe 


Z 


11 


18 


32 


53 


54 


1*1/2 


-40.54 


-119.1 


-411.0 


-1225. 


-1277. 


2Sl/2 


-2.805 


-12.41 


-53.45 


-193.0 


-202.4 


2pi/2 


-1.522 


-9.631 


-47.33 


-180.5 


-189.6 


2^3/2 


-1.514 


-9.547 


-46.14 


-169.5 


-177.7 


3S1/2 


-0.1823 


-1.286 


-7.409 


-40.52 


-43.01 


3pi/2 




-0.5953 


-5.324 


-35.34 


-37.66 


3^3/2 




-0.5878 


-5.157 


-33.21 


-35.32 


3*^3/2 






-1.616 


-24.19 


-26.02 


3^5/2 






-1.591 


-23.75 


-25.53 


4si/ 2 






-0.5687 


-7.759 


-8.430 


4P1/2 






-0.2821 


-5.868 


-6.452 


4P3/2 






-0.2730 


-5.450 


-5.982 


4d3/2 








-2.341 


-2.711 


4d 5 /2 








-2.274 


-2.633 


5S1/2 








-0.8762 


-1.010 


5P1/2 








-0.4341 


-0.4925 


5P3/2 








-0.3903 


-0.4398 



Kis = fi(Z, e a + e ls ) 
K 2s = f2(Z,e a + e 2s ) 



3847^1, 



-4ui arccoti/i 



(e Ze a ) 2 l- e -2^i 



61447re| 



1 + 3^ 



K 2p = 



g— Av2 arccot (1^2/2) 
X n , 

12288^e 4 

J 2 (Z, e a + e 2p )- 



3 + 8— 



-4^3 arccot(z^3 /2) 

1 _ e -27TI/ 3 



(10) 



(11) 



(12) 



where a is the fine structure constant, Z is nuclear 
charge, e a is axion en ergy, e 2 = e 2s |, e 3 = |e 2p |, 
V\ = \f-e ls /(e ls + e a ), v 2 = 2^/-e 2s /(e 2s + e a ), z/ 3 = 
2^/—e 2p /(e 2p + e a ). Here ei s , £2 S and e 2p are the Hartree- 
Fock energies of the core states. Hartree-Fock energies of 
all core states of Na, Ar, Ge, I and Xe are presented in 
Table Q] For other atoms extrapolation formulas can be 
used: 



(Z) = - 



en 
Q) 
^i(Z) 



^E(Z) 



- 7.49Z + 43.39 

2 ' 

-0.000753Z 3 - 0.028306Z 2 
-0.066954Z +2.359052, 
-0.000739Z 3 - 0.027996Z 2 
+0.128526Z + 1.435129. 



(13) 

(14) 
(15) 



The functions fi(Z) and f 2 (Z) in (|1QI11I12[) are scaling 
functions: 



fx(Z,e) 



(5.368 x IQ- J Z - 1.17 X 10" 4 )e/e 
-0.012Z+ 1.598 (16) 
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FIG. 1: Dimensionless factor K (see formula ((3| in the ion- 
ization cross sections of Ar by axion. Solid line - result of 
Hartree-Fock calculations, dotted line - formula ((9|. 



f 2 (Z,e) = (-1.33 x 10~ 6 Z + 1.17 x 10" 4 )e/e 

-0.0156Z+ 1.15 (17) 

They are chosen to fit the results of Hartree-Fock cal- 
culations for axion energies between 10 and 50 keV and 
nuclear charge 18 < Z < 60. We stress once more that 
not only scaling functions but all the formulae (|10H17j) 
were obtained by fitting the Hartree-Fock calculations. 
All atomic shells are included into the RHF calculations. 
Therefore, the fit with K to tai = K\ s +K 2s + K 2p includes 
the higher shells as well. As it will be demonstrated in 
the next section, the analytical and Hartree-Fock results 
agree within few per cent. 



III. RESULTS 



Figures II 121.31 show the results of relativistic Hartree- 
Fock and analytical calculations of the dimensionless fac- 
tor K which stands in the expression for the cross section 
of atom ionization by axion (see formula ©). We use 
the RHF computer code and formulae ^ and (j3]) for the 
calculations. Many body and relativistic effects beyond 
the RHF method are ignored and the final electron state 
in the continuum is calculated in the same potential as 
initial core state. The accuracy of this approximation 
is few percents due to dominating contribution from the 
inner-most core states Is, 2s and 2p. For these states the 
many-body effects are small due to strong nuclear field. 

Calculations using analytical formula (|9I10I11I12[) are 
shown on Figures I1I2I3I as dotted lines. The agreement 
between RHF and analytical results is very good for all 
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FIG. 2: As on Fig. [Lj but for Ge. 
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three atoms. It is on the order of a few per cent when an 
electron from the Is is excited and better than 20% when 
only 2s and higher states contribute to the cross section. 
We expect similar trend for other atoms. This means 
that at present stage the analytical formulae (j9H0ll 1I12[) 
are sufficiently accurate and there is no need for more 
sophisticated atomic calculations. 

The results show an obvious advantage of using heavy 
elements in the detectors given that axion energy is suffi- 
ciently large to excite an electron from an inner-most sub- 
shell. For example, for a solar axion of e a =14.4 keV [19j . 
K — 450 for the case of Ar and K m 4000 for the case 
of Ge and Xe. Current experimental sensitivity to the f a 
parameter is in the range of 10 s — 10 9 GeV 0, Q, but is 
poised to be improved in the very near future. 

In this paper we have examined the case of atomic 
ionization by the pseudoscalar super-weakly interacting 
dark matter, but a similar case can be made for other 
types of super- WIMPs with arbitrary integer spin, and 
coupled to electrons via vector, axial-vector, tensor etc 
couplings. Another obvious use of the RHF calculations 
relevant to the underground searches of exotic particles, 
is the generalization to the detection of nearly massless 
particles emitted by the solar interrior with typical ener- 
gies of a few keV. Both applications will be considered in 
a forthcoming publication. 
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FIG. 3: As on Fig. CD but for Xe. 
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